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Photonic crystals doped with resonant atoms allow for uniquely advantageous nonlinear modes of 
5_l , optical propagation: (a) Self-induced transparency (SIT) solitons and multi-dimensional localized 

Oh "bullets" propagating at photonic band gap frequencies. These modes can exist even at ultraweak 

intensities (few photons) and therefore differ substantially either from solitons in Kerr-nonlinear 
photonic crystals or from SIT solitons in uniform media, (b) Cross-coupling between pulses exhibit- 
ing electromagnetically induced transparency (EIT) and SIT gap solitons. We show that extremely 
strong correlations (giant cross-phase modulation) can be formed between the two pulses. These 
features may find applications in high-fidelity classical and quantum optical communications. 
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Photonic crystals (PCs) can exhibit an interplay between Bragg reflections, which block the propagation of light in 
photonic band gaps (PBGs ) ffl [|, |6), and the dynamical modifications of these reflections by nonlinear light- 

matter interactions ]j, ^, |10|, |ll|, [12["l A very interesting situation arises when foreign atoms or ions — dopants — with 
O |' transition frequencies within the PBG are implanted in the PC ||, Then light near one of these frequencies 
resonantly interacts with the dopants and is concurrently affected by the PBG dispersion. Consequently, highly 
£h nonlinear processes with a rich variety of unusual PC-related features are anticipated. 

Our aim in recent years has been to identify those regimes of nonlinear optical propagation in doped PCs that allow 
^ transmission of extremely weak pulses, while filtering out undesirable noise, and are therefore highly advantageous for 
optical communications, data storage and processing, near or at the quantum limit. These requirements are satisfied 
by novel regimes surveyed in this paper that have been theoretically discovered and investigated by us: a) Self-induced 
transparency (SIT) solitons propagating inside or near a PBG at a frequency that is near resonant with the transition 
frequency of the dopant jl3|, [l4|, [l5| [l(|. This peculiar form of gap solitons (GSs) is immune to resonant absorption 
even for a small number of photons, and may also possess two- or three-dimensional (2D or 3D) localization in the 
form of light bullets (LBs)]l7| (Sec. II). b) Cross- coupling of SIT and electromagnetically-induced transparency (EIT) 



pulses in PCs. We put forward a new regime in Sec. [II: a strong modulation of the phase of a weak pulse subject 
to EIT by a control pulse in the form of an SIT GS moving at the same slow velocity. Thereby giant cross-phase 
modulation can be formed between the GS and the EIT pulses. In Sec. |y| we summarize and discuss our findings in 
this paper. 



II. SELF-INDUCED TRANSPARENCY (SIT) GAP SOLITONS 



A. Background 



A GS is usually understood as a self-localized moving or standing (quiescent) bright region, where light is confined 
by Bragg reflections against a dark background. The soliton spectrum is tuned away from the Bragg resonance 
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by the nonlinearity at sufficiently high field intensities. The first type of GS had been predicted |18|, [jjj [2C], 
and later observed |p2| , in a Bragg grating possessing Kerr-nonlinearity. A principally different mechanism of GS 
formation has been theoretically discovered by our group in a periodic array of thin layers of resonant two-level 
atoms (TLA) separated by half- wavelength nonabsorbing dielectric layers, i.e., a resonantly absorbing Bragg reflector 




15| , As opposed to the 27r-solitons arising in SIT, i.e., resonant field-TLA interaction in uniform 
their GS counterparts in a RABR may have an arbitrary pulse area JIJ, 15 1. It must be stressed 



(RABR 
media 

that stable, moving or standing, GS solutions have been consistently obtained only in a RABR with thin active TLA 
layers. By contrast, the case of a periodic structure uniformly doped with active TLA calls either for a solution of the 
wave (Maxwell) equation without the spatial slowly- varying envelope approximation (SVEA), or for a solution of an 
infinite set of coupled Bloch equations for all spatial harmonics of the atomic polarization (Fourier components) [ jl3| . 
Therefore, an attempt [ p5| to obtain a self-consistent solution for a uniformly-doped periodic structure by imposing 
the SVEA, or by arbitrarily truncating the infinite hierarchy of equations for the harmonics of atomic population 
inversion and polarization to its first two orders, is generally unjustified. In fact, it can be shown numerically to fail 
for many parameter values. 

In the simplest case of a uniform (bulk) medium, when the driving field is resonant with the atomic transition, the 
TLA Bloch equations can be easily integrated and the Maxwell equation then reduces to the sine-Gordon equation 



d 2 e . 
dCdf = - sm0 

for the pulse area 9 — f_ Qdt', i.e., the time integral of the Rabi frequency f2. Equation 
the dimensionlcss variables f — (t — noz/c)/To and £ = hqz/ctq, where tq = ^ 



(1) 

is written in terms of 
is the cooperative resonant 

absorption time, go being the TLA density (averaged over z), fi the dipole moment of the TLA transition at frequency 
wo and no is the refraction index of the host medium. This sine-Gordon equation is known to have solitary-wave 
solutions, which propagate without attenuation or distortion with a conserved pulse area of 2ir |p3| , 24 1. These SIT 
solitons have the form: 



n((,f) = (T )- 1 A sech[f3((-vf)}, 



(2) 



where the pulse width (3 is an arbitrary real parameter uniquely defining the amplitude Aq — 2//3 and group velocity 
v = 1//3 2 of the soliton. In what follows, Eq. (||) will be compared with an SIT GS in a RABR. 



B. SIT in RABR: The Model 

Let us assume |l3|, O, 15, [Ig| a one-dimensional (ID) periodic modulation of the linear refractive index n 2 (z) = 
7ig[l + a\ cos(2fc c z)]. The periodic grating gives rise to a PBG with a central frequency uj c — k c c/no and gap edges 
at wi,2 = 07 c (l±ai/4). The electric field £ of a pulse propagating along z can be expressed by means of the 
dimensionless quantities T,± = 2r /z?i _1 (£p ±£b), where £f and £b denote the slowly varying amplitudes of the 
forward and backward propagating fields, respectively, as 

E(z,t) = h(fiTo)- 1 {Re[S + (z,i)e~ lw<:t ]cosfc c z - Im[£_(z, i)^^*] sin£; c z} . (3) 

We further assume that very thin TLA layers (much thinner than l/k c ), whose resonance frequency ujq is close to the 
gap center oj c , are placed at the maxima of the modulated refraction index (Fig. |l|). They are located at positions Zj 
such that 



e 



ik c Z2j i ~ik 



1, e 1 ^ 2 ^ 1 = -1, (4) 



i.e., the TLA density is described by g — gaX/Z^j 8{z — Zj), where A = 2Tr/k c is the wavelength. 

The Bloch equations for the slowly varying polarization envelope P and inversion w in the even numbered 
layers can be obtained (in the slowly varying envelope approximation) by substituting for the Rabi frequency 
£1 = Tq 1 ('E + cos k c z + iS_ sinfc c z) and applying Eq. (Q) at the positions of these layers: 

dP 

— = -iSP + Z+w, (5) 
^ = -Re (£ + P*). (6) 
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FIG. 1: Schematic description of the periodic RABR and of the decomposition of the electric field E into modes £+ and £_. 
The shading represents regions with different index of refraction; the darker the shading the larger n is. The black regions 
correspond to the TLA layers. The localization of the field envelope over ~ 20 structure periods is shown for the sake of 
visualization; in reality, the field is localized over a hundred or more periods. 



Combining Eqs. (|^) and (^), one can eliminate the TLA population inversion: w — s/i - \P\ 2 . The remaining 
equation, together with the Maxwell equations for £± (driven by P), form a closed system, 

7 1 2 Z + + 2i(r 1 -5)P-2^/l-\P\^ + , (7) 
BP 

VS--2|, (8) 

-iSP- v/l-|P| 2 S+, (9) 

where r = £/tq, £ = (uq/ctq) z and 6 = (ujq — w c )to are the dimensionless time, coordinate, and detuning, respectively, 
and T) = /abs/^rcfl = aiLu c T /4: is the dimensionless modulation strength, which can be expressed as the ratio of the 
TLA absorption distance l^bs — tqc/uq to the Bragg reflection distance / re fl = 4c/(aiw c no). We emphasize that the 
above equations are obtained using the SVEA, which is valid under the assumption that the Bragg reflection does 
not appreciably change the pulse envelope over a distance of a wavelength, i reH 3> A, whence a% <C 2/7r. 

To reach general understanding of the dynamics of the model, one should first consider the spectrum produced by 
the linearized version of Eqs. which describes weak fields in the limit of infinitely thin TLA layers. Setting 

£+ = Ae l{K ^ XT \ £_ = Be'(^^"ffi = — 1, and P — Ce l{ -^-* T \ we obtain from the linearized equation (§) that 
C = i(6 — x) ~ 1 A. Substituting this into Eqs. (0) and (|^), we arrive at the dispersion relation for the wavenumber k 
and frequency x m the form 

(X 2 - - rf)( X -S)x{( x -S) [ X 2 - n 2 - (2 + V 2 )] + 2(r? - 6)} = 0. (10) 

Different branches of the dispersion relation generated by Eq. ( |l0| ) are shown in Fig. ^. The roots x — i\/ k 2 + rj 2 
(corresponding to the solid lines in Fig. ^) originate from the driven equation (||) and represent the dispersion relation 
of a Bragg reflector with the gap |x| < Vi that does not feel the interaction with the active layers. Important roots 
of Eq. ( |l0|) are those of the expression in the curly brackets, shown by the dashed and dash-dotted lines in Fig. [I]. 
These roots correspond to nontrivial spectral features: bright or dark solitons in the indicated (shaded) bands. 

The frequencies corresponding to k = are Xo = V and Xo,± — ~ (v ~ S)/2 ± \/2 + (rj + 5) 2 /2, while at k 2 — > oo 
the asymptotic expressions for different branches of the dispersion relation are x — i& an( i X — $ + 2 (rj — S) k~ 2 . 
Thus, the linearized spectrum always splits into two gaps, separated by an allowed band, except for the special case, 
X] = ?7o = 6/2 + \Jl + S 2 /i, when the upper gap closes down. The upper and lower band edges are those of the 
periodic structure, shifted by the induced TLA polarization in the limit of a strong reflection. They approach the SIT 
spectral gap for forward- and backward-propagating waves (2(| in the limit of weak reflection. The allowed middle 
band corresponds to a polaritonic excitation (collective atomic polarization) in the periodic structure. 
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FIG. 2: The RABR dispersion curves at r\ = 0.5 and S — —0.2. The solid lines show the dispersion branches corresponding 
to the 'bare' (noninteracting) grating, while the dashed and dash-dotted lines stand for the dispersion branches of the grating 
'dressed' by the active medium. The frequency bands that support the standing dark and bright solitons are shaded. 




C. Standing (quiescent) self-localized pulses 

We seek the stationary solutions of Eqs. (Q) and ^ corresponding to bright solitons in the form 



with real V and S. Substituting this into 



= e-' ixr «S(C), P = ie- lXT T(Q (11) 
3), we eliminate V in favor of S and obtain an equation for <S(£), 
(V ~ X) ■ sign(x - 5) 



(12) 



It then follows |L5j that bright solitons can appear in two frequency bands \t the lower band being xi < X < 
min{x2, — f], S}, and the upper band being max{xi, r), 5} < x < X2, where \i,2 = 1/2 [(5 — r\ =F \J(j] + &) 2 + 8] are 
the boundary frequencies. The lower band exists for all values rj > and 8, while the upper one only exists for the 
weak-reflectivity case S > i] — 1/rj. On comparing these expressions with the spectrum shown in Fig. pi we conclude 
that part of the lower gap is always empty from solitons, while the upper gap is completely filled with stationary 
solitons in the weak-reflectivity case, and completely empty in the opposite limit. 
In an implicit form, the solution of Eq. (113) reads 



s(0 = 2\x-s\n(0(i-n 2 (0) 1 



with 



ICI 



x 



[l-TZly 1/2 tan 



IK* 



TV 



i-nl 



(2^ )- 1 ln 



K 



(13) 



(14) 



and TZq = 1 — | (x + 77) (x — (5) | /2. This zero- velocity (ZV) gap soliton is always single-humped and its amplitude, found 
from Eq. (0), is given by 



^o/Vlx + vl 



(15) 



To calculate the electric field in the antisymmetric £_ mode, we substitute £_ = ier 1 * 8 A{Q into Eq. (|) and 
obtain 



A" + (x 2 ~ V 2 ) A = 2V', 



(16) 



which can be easily solved by the Fourier transform, once V(C) is known. We note that, depending on the parameters 
rj, 6 and %, the main part of the soliton energy can be carried either by the S + or the S_ mode. 




FIG. 3: Pulses obtained as a result of 'pushing' a zero- velocity RABR soliton (dashed lines): (a) push, characterized by the 
initial multiplier exp(— ip() after a sufficiently long evolution (r = 400) (solid lines). 5 = 0, n — 4, \ — —4.4, and p = 0.1. (b) 
idem, but for p — 0.5. 



The most drastic difference of these new solitons from the well-known SIT solitons in Eq. (g) is that the area of 
the ZV soliton (integrated over £) is not restricted to 2tt, but, instead, may take an arbitrary value. This basic new 
feature shows that the Bragg reflector can enhance (by multiple reflections) the field coupling to the TLA, so as to 
make the pulse area effectively equivalent to 2ir. In the limit of the small-amplitude and small-area solitons, TZq -C 1, 
Eq. ( [j~4] ) can be easily inverted, the ZV soliton becoming a broad sech-like pulse: 



n (j . (17) 

In the opposite limit, 1 — TZq — > 0, i.e., for vanishingly small \\ + f]\, the the soliton is characterized by a broad central 

part with a width ~ (l — TZq) and its amplitude ( |l5| ) becomes very large. Thus, although the ZV soliton has a 
single hump, its shape is, in general, strongly different from that of the traditional nonlincar-Schrodinger (NLS) sech 
pulse. 



S^2\ x -S\TZq sech 



x-v 



D. Moving solitons 

One could expect a translational invariance of the ZV solitons (|l|) on physical grounds. Hence, a full family of 
soliton solutions should have velocity as one of its parameters. This can be explicitly demonstrated in the limit 
of the small-amplitude large- width solitons [cf. Eq. Jl7|)]. We search for the corresponding solutions in the form 
£+(C> t ) = ^(Cj t ) ex P( — iXo T ), P (C; r ) — ^(C) r ) ex P( — *Xo T ) [cf- Eqs. (|il])1, where xo is the frequency corresponding 
to k — on any of the three branches of the dispersion relation ( [To| ) (see Fig. and the functions S((, t) and V{C, 7 r) 
are assumed to be slowly varying in comparison with exp(— ixo T )- Under these assumptions, we arrive at the following 
asymptotic equation for S{Q, r): 

5 =(" 2 -*" + 2 S^) 5 - (18 » 

Since this equation is of the NLS form, it has the full two-parameter family of soliton solutions, including the moving 
ones 0- 

In order to check the existence and stability of the moving solitons numerically, the following procedure has been 
used |ll|: Eqs. (Q) and ([)]) were simulated for an initial configuration in the form of the ZV soliton multiplied by 
exp(ipC) with some wavenumber p, in order to 'push' the soliton. The results demonstrate that, at sufficiently small 
p, the 'push' indeed produces a moving stable soliton [Fig. ||(a)]. However, if p is large enough, the multiplication by 
exp(ipC) turns out to be a more violent perturbation, splitting the initial pulse into two solitons, one quiescent and 
one moving [Fig. 13(b)] . 
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FIG. 4: The forward-propagating electric field of the two-dimensional 'light bullet' in the Bragg reflector, \£f\, vs. time r and 
transverse coordinate x, after having propagated the distance z = 1000. The parameters are n = 0.1, S = 0.2, C = 0.1 and 
Oo = —1000. The field is scaled by the constant h/4rofj,no. 



E. Light bullets (spatiotemporal solitons) in PCs 



The advantageous properties of SIT GS can be supplemented by immunity to transverse diffraction, i.e., simultane- 
ous transverse and longitudinal self-localization of light in a PC: multi-dimensional spatio-temporal solitons or "light 
bullets" (LBs) |28| have been analytically and numerically predicted by our group to exist and be stable, not only in 
uniform 2D and 3D SIT media ^9|, but also in 2D or 3D periodic structures, wherein SIT solutions combining LB 
and GS properties are demonstrated [ fl7| . Our objective is to consider the propagation of an electromagnetic wave 
with a frequency close to lu c through a 2D PC doped by thin TLA layers. The forward- and backward-propagating 
components satisfy equations that are a straightforward generalization of the ID equations (0) and (||) 

. d 3 E + | ; d 3 £- | d 2 £+ d 2 £+ | ^ d 2 £+ | 2£ _ 2 ®L-2iP~0 (19) 
drx 2 d^x 2 dr 2 dQ 2 dx 2 dr 
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h dr 2 



. d 3 E„ | , d 3 E + | d 2 E, d 2 T,^ ^ d 2 E_ | 2 ^ + 2 dP_ Q , 2Q , 
drx 2 d^x 2 dr 2 dC, 2 dx 2 dQ 

where the Fresnel number F > 0, which governs the transverse diffraction in the 2D and 3D propagation, was 
incorporated into x denoting the transverse coordinate. The equations for the polarization P and inversion w are the 
same as Eqs. (||) and 

We search for analytical LB solutions of Eqs. (jl9|), (po|), (||) and (||), by the following ansatz that reduces in ID to 
the exact moving GS O, 111 



£+ = A Vsecheiseche 2 e 4(KC - XT)+ " /4 , (21) 
= £+/«, (22) 
P = v/secheiseche 2 {(tanhei + tanh9 2 ) 2 + 

^^C 4 [(tanh 9i - tanh 9 2 ) 2 - 2(sech 2 9i + sech 2 2 )l 2 }V V^-*^"', ( 2 3) 
4rj 

- = [l-l^l 2 ] 172 , (24) 

with 9i(r, () = (](( — vt) + 9o + Cx, 6 2 (r, C) = P{( — vt) + 9o — Cx, the phase v and coefficients 9o and C being 
real constants, while the other parameters are denned as A = 2^/6/^-1, [3 = y/S/r]+l, v = - y/{8 - r])/(5 + r/), 
K = —\/5 2 — r] 2 , and % = <5. 

The ansatz (21)- ( p4| ) satisfies Eqs. ( pj| ) and ( pp| ) exactly, while Eqs. (Q) and (||) are satisfied to order ^S/rj — IC 2 , 
which requires that y/6/r] — 1C 2 <C 1. The ansatz applies for arbitrary 77, admitting both weak (77 <C 1) and strong 
(r; > 1) reflectivities of the Bragg grating, provided that the detuning remains small with respect to the gap frequency. 
Comparison with numerical simulations of Eqs. (|l9|), (|20|), (||) and (g), using Eqs. (^l])-(|24|) as an initial configuration, 
tests this analytical approximation and shows that it is indeed fairly close to a numerically exact solution; in particular, 
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FIG. 5: (a) Schematic representation of the atomic system interaction with strong driving field on the transition |2) — » |3) and 
weak fields E a and on the transitions |1) — > |2) and |3) —* 4), respectively, (b) Absorption and dispersion spectra of the 
E a field in the absence (solid line) or presence (dashed line) of the Eh field. 



the shape of the bullet remains within 98% of its originally presumed shape after having propagated a large distance, 
as is shown in Fig. |[ 

Three-dimensional (3D) LB solutions with axial symmetry have also been constructed in an approximate analytical 
form and successfully tested in direct simulations, following a similar approach ]l7| . Generally, they are not drastically 
different from their 2D counterparts described above. 



F. Information transmission by SIT GSs and LBs 



The efficiency of information transmission is characterized either by channel (information) capacity C = 
Wln(I s /I n ), where W is the bandwidth and I s /I n is the ratio of the signal-to-noise intensities, or by the data 
transmission density D = NM, where N is the number of bits per channel and M is the number of accessible chan- 
nels. Both C and D can be very high in the case of a SIT GS or LB for the following reasons: (a) The bandwidth 
W is large, being limited by the PBG width, which can be very large in the optical domain. At the same time, noise 
is very effectively suppressed by the Bragg reflection and by the absence of diffraction losses in the case of a LB. (b) 
The maximal transmission density D can be estimated |30| as the ratio of the accessible bandwidth, in our case the 
PBG width (in excess of 10 13 s _1 in the optical domain), to the spontaneous linewidth (10 6 s _1 for rare-earth ions). 
Hence, these modes of transmission can be very effective for optical communications. 



III. CROSS COUPLING BETWEEN ELECTROMAGNETIC ALLY-INDUCED AND SELF-INDUCED 

TRANSPARENCY PULSES 



A. EIT in bulk media: Background 



Electromagnetically induced transparency (EIT) is based on the phenomenon of coherent population trapping 
|32f , in which the application of two laser fields to a three-level atomic system creates the so-called "dark state" , 
which is stable against absorption of both fields. Consider a four- level atomic system whose level configuration is 
depicted in Fig. |5|(a). The phase shift and absorption of an optical field Ei are given by the real and imaginary parts 
of its complex polarizability a«. In the absence of the "control" field Et, the usual EIT spectrum [Fig. |^(b)] for the 
weak probe field E a exhibits vanishing phase shift and absorption [Re(a a ) = Im(a a ) — 0] at the two-photon Raman 
resonance oj a = ui^x + uJd, where uo a and Ud are the frequencies of the probe and driving fields, respectively, and uiij 
is the frequency of the atomic transition \i) — > \j). An off-resonant control field E^ with the frequency u>b such that 
\At \ = \b)b — 3> 74, where 7* is the decay rate of the corresponding atomic level, induces an ac Stark shift of level 
|3) and thereby shifts the EIT spectrum [Fig. |](b)]. 

Due to the steepness of the dispersion curve in the vicinity of the Raman resonance, |(9 Wa Re(a! a )| ^ |9 Wa Im(a(,)|, 
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ID |3> 15) 

FIG. 6: (a) Schematic representation of the field-atom system: The probe field E a exhibits EIT on the transition |1) — * |2), 
in the presence of the driving field Ed on the transition |3) — > |2); The control field Eb exhibits SIT GS on the transitions 
5) — > |6); The transition |3) — > |4) serves to cross-couple the two fields. Initially only the states |1) and |5) are populated, 
(b) Polarizations and propagation directions of the fields involved in the transitions above. 



this Stark shift leads to a large phase shift along with small absorption of the probe field E a : 

, x> f \ ao72|^b| 2 j , , 7 4 Re(a a ) 

<f> a = Re{a a )z ~ - — z, Im(a a ) = — < Re(a a ), (25) 

where ao is the resonant absorption coefficient of the medium at the frequency W21 and Oj = /i^Ei/h is the Rabi 
frequency of the corresponding field (/ijj the dipole matrix element on the respective transition). This is the essence 
of the so-called giant Kerr cross-phase modulation of a probe field by a control field, introduced first by Schmidt and 
Imamoglu [Q. Later Harris and Yamamoto |p3 have predicted that a resonant control field Eh with \Ah\ < 74 can 
destroy the coherence between the two ground levels |1) and |3), which leads to a two-photon absorption Im(a a: b) = 

i2 

that is, the medium absorbs two fields simultaneously, but does not absorb one field alone. This is the 



ao72 \£h 



essence of a probe-photon switch, conditional on the presence of control photons. 

The main limitation of the above schemes |33|, |34|, |3^] stems for the fact that the effective interaction length is 

limited by the mismatch between the group velocity of the slowly propagating E a field, ~ <C c/n and that 

of the nearly-free propagating Eh field, Vg ~ c/uq. For weak (few-photon) pulses, this mismatch ultimately limits 
the maximal phase shift or absorption of the probe in the presence of the control field. 



B. Simultaneous EIT and SIT in RABR 



In this section we propose a new implementation of the cross-phase modulation, in which the group velocities of 
both fields can be matched, allowing one to obtain any desired phase shift of the probe field with a weak control field. 
To this end, we consider the same configuration as in Sec. II, leading to SIT GS and LB solutions: a PC periodically 
doped by thin layers of atoms at the maxima of its refractive index. However, the multi-level structure of the atoms 
is now playing a role: it is shown in Fig. ^|, along with the polarizations and propagation directions of the fields 
involved. The states |1), |3) and |5) are the degenerate Zeeman components with Mp = — 1,0, +1, respectively, of 
the atomic ground level having total angular momentum F = 1. Similarly, the states |4) and |6) are the degenerate 
Zeeman components with Mp = —1,0, respectively, of the excited level having angular momentum F = 1. Finally, 
the state \2) corresponds to the single Zeeman component with Mp = of another excited level having F = 0. Such 
a level scheme is found, e.g., in alkali atoms, where the ground level is Si/%,F = 1 and the two excited levels are 
P1/21F — 1 an d Pz/2,F = 0. Due to the dipole selection rules, the 7r-polarized driving field couples the states with 
AM = 0, the cr+-polarized E a field couples the states with AM = 1 and the er_-polarized Eh field couples the states 
with AM = -1. 

We assume that initially all the atoms are optically pumped into the states |1) and |5), which then acquire equal 
populations 1/2. Hence, the sequence of trans itions |1) — > \2) — > |3) — » |4) repeats that of Fig. pRa), realizing the 



cross-phase modulation scheme of Sec. Ill [II A. The frequency of the E a field is far from the band gap frequencies of 



the PC, while the frequency of the Eb field is within the band gap. 
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FIG. 7: Schematic representation of the level scheme suitable for EIT in the presence of a Raman GS (only the relevant levels 
are shown) 



As was shown in Sec. II, PC structures doped with the near-resonant TLAs can support standing and slowly 
moving SIT GSs, whose pulse area (integrated over z) can take an arbitrarily small value. In the present setup, the 
transition |5) — > |6) realizes that near-resonant TLA, allowing for slow propagation of the Eb field through the PC. 

Let us write the propagation equation for the slowly moving SIT soliton in the form 

(tv {h) -z\ 

fl slT (z, t) = ft b sech I 9 j , (26) 

where fib is the peak Rabi frequency, v g is the group velocity and (3 = (2ctt>) _1 , with at, being the absorption 
coefficient of the active medium at the carrier frequency u>b of the soliton. The temporal width of the pulse is given 

(b 



by Tb = 2/3/vg b \ The area of the Eb pulse 



Q slT dt = -^-tt, (27) 

Vg OLb 



is then inversely proportional to the group velocity of the pulse. Hence, the SIT condition (Sec. [I HII A ) 9b = 2-k 
imposes a unique relation between the Rabi frequency of the SIT soliton and its group velocity: 



(b) _ 
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h 



2a b 



(28) 



The absorption- free propagation of the SIT soliton is limited to z < /T6; where is the decay rate of the upper 
atomic state |6). 

Our aim is to match the group velocities of the E a field subject to EIT and the Eb field having the form of a slow 
SIT gap soliton: Vg = v g b \ This requires that \fld\ 2 = flbCt a l2/ ^otb, i.e., an appropriate choice of the driving field 
Rabi frequency fid, for a given Rabi frequency fib of the control field. Such velocity matching of the two copropagating 
weak fields would maximize their interaction. 

One possibility to launch the required slow SIT soliton is to irradiate the PC by a laser beam at a small angle 
tp relative to the periodicity direction z, ip — DrioVg^ / (Lc) <C 1, where D and L are, respectively, the transverse 
and longitudinal dimensions of the structure. This choice of ip ensures that the z-component of the beam, which 
forms the SIT soliton and propagates in the PC with the group velocity v g over the distance L, will traverse the 
structure during the same time as the transverse component of that beam, which covers the distance D with the 
velocity (c/no) sin?/\ 

We have checked that for the parameter values corresponding to dopant atoms (or ions) with the mean density 
N = 10 13 cm" 3 (surface density of 4 x 10 8 cm" 2 in the thin layers), A fc = 3O74, |O b | ~ 10 6 rad/s and \ fl d \ ~ 4 x 10 6 
rad/s, we obtain tt phase shift of the E a field over a distance z ~ 4 cm, while the absorption probability remains less 
that 10%. 

One possible difficulty of our scheme is that, with the parameters listed above, the temporal width of the Eb field 
is Tb ~ 10~ 6 s, and the interaction time is of the order of 10 -3 s, while the lifetime of the SIT soliton is of the order 
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of the decay time of the excited atomic state 7^ ~ 10~ 7 s. One can cope with this problem by employing the atomic 
level scheme shown in Fig. [z| which allows one to launch Raman solitons. We irradiate the system with an additional 
strong cw field Er, which couples the F = 1 excited state with the F = 2 metastable ground states. The fields Er 
and Ei, are largely detuned from the fast decaying excited states |4) and |6) by an amount A 3> 74,6- Then, upon 
adiabatically eliminating the states |4) and |6), we obtain that the Rabi frequency fib of the control field in Eq. ( prjj ) 
is simply replaced by fibOfl/A. The lifetime of the SIT soliton is given now by the lifetime of the F — 2 ground 
states, which can be very large, reaching in some instances a fraction of a second! In addition, such a setup allows 
one to launch slow Raman GSs []36[ , and thus circumvent the difficulty of launching a standing (ZV) or slowly moving 
GSs, which must otherwise overcome the high reflectivity of the PC boundaries. 



IV. CONCLUSIONS 



In this paper we have focused on properties of solitons in a doped PC or RABR, combining a periodic refractive- 
index superlattice (Bragg reflector in ID or 2D) and a periodic set of thin active layers (consisting of TLAs resonantly 
interacting with the field). We have demonstrated that the system supports a vast family of bright GSs, whose prop- 
erties differ substantially from their counterparts in periodic structures with either cubic or quadratic off-resonant 
nonlinearities. Depending on the initial conditions, these can be either standing (ZV) or slowly moving stable solitons 
that exhibit SIT irrespective of their photon number (pulse energy) for an appropriate group velocity. A multidi- 
mensional version of this model corresponds to a periodic set of thin active layers placed at the maxima of a 2D- or 
3D-periodic refractive index. It has been found to support stable propagation of spatiotemporal solitons in the form 
of 2D- and 3D-localized LBs. 

The best prospect of realizing a PC which is adequate for observing the GSs and LBs is to use thin layers of 
rare-earth ions |37j embedded in a spatially-periodic semiconductor structure p8[ . The TLAs in the layers should be 
rare-earth-ions with the density of 10 15 — 10 16 cm -3 , and large transition dipole moments. The parameter rj can vary 
from to 100 and the detuning is ~ 10 12 — 10 13 s" 1 . Cryogenic conditions in such structures can strongly extend the 
dephasing time T2 and thus the soliton's or LB's lifetime, well into the fisec range |37j, which would greatly facilitate 
the experiment. 

In a 2D PC, LBs can be envisaged to be localized on the time and transverse-length scales, respectively, ~ 10~ 12 
s and 1 /im. The incident pulse has uniform transverse intensity and the transverse diffraction is strong enough. 
One needs d 2 // a b s Ao < 1, where ^ a b s , Ao and d are the resonant-absorption length, carrier wavelength, and the pulse 
diameter, respectively. For Z a b s ~ 10~ 3 m and Ao ~ 10~ 4 m, one thus requires d < 10~ 4 m, which implies that the 
transverse size of the PC must be a few /i m. 

We have considered here (Sec. Ill [II B) the cross-coupling of optical beams in a PC or RABR. We have pointed 
out, for the first time, the advantageous features of the cross coupling between EIT and SIT pulses, which is capable 
of producing extremely strong correlations between the two pulses. With doping parameters as above, and driving 
and control fields with Rabi frequencies of the order of 10 6 rad/s, we can obtain a phase shift of it for the weak probe 
pulse over a distance of few cm. This is much larger than any corresponding phase shift (for similar control fields) in 
other media. 

We strongly believe that the highly promising payoff expected from the construction of suitable structures justifies 
the experimental challenge they pose. If and when the schemes proposed above are experimentally realized, they may 
prove to be useful for producing ultrasensitive nonlinear phase shifters or logical photon switches for both classical 
and quantum information processing or communication, owing to the unique advantages of the doped PCs over 
conventional EIT schemes §| |§ || or high-Q cavities @ 
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